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ABSTRACT 

Quadrotor control is needed so that the quadrotor ean float close to the stationary 
state. For that we need control techniques. One control technique that can be 
designed and implemented in quadrotor is PID eontrol. PID parameter tuning using 
the Genetic Algorithm technique can speed up the manual tuning process. The 
weakness in the applieation of the Genetic Algorithm rule is that it often rejeets 
important information found in other individuals and causes premature convergence, 
especially at the beginning of the generation. These problems can be overeome by 
using crossover and mutation rules with different probability levels according to 
fitness values and evolutionary processes. The results of the study using fast genetie 
algorithm techniques obtained constants Kp, Ki and Kd with the lowest rise time and 
overshoot, namely 0.010, 0.001 and 0.036 at the pitch angle. At the roll angle, they 
are 0.010, 0.001 and 0.03. At yaw angle 0.018, 0.006 and 0.043. Comparison of PID 
tuning simulations using fast genetie algorithm with genetic algorithm standards, 
shows that fast genetic algorithm has increased optimum generation achievement 
faster by 26.67% at pitch angle, 44%) at roll angle and 20% at yaw angle. This 
condition has an effect on increasing simulation execution time, where fast genetic 
algorithm is 26.4% faster at piteh angle, 38.05% at roll angle, and 24.19% at yaw 
angle. 
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INTRODUCTION 

Quadrotor or often called quadcopter is a type of Unmanned Aerial Vehicle (UAV). 
Quadrotor can fly in the air using the help of four motors and four propellers (Dharmawan 
and Putera, 2012). By using the thrust of the rotation of each motor, it can produce movement 
on the quadrotor. The control system that regulates this is the attitude control system. The 
attitude control system can be generated through system modeling. At present there are two 
modeling methods that can be used for attitude from quadrotor dynamics, namely Lagrange 
equations and Newton II laws (Dharmawan et al., 2014). As for the well-known control 
system and can be implemented directly on quadrotor, it is a PID control (Bouabdallah et al., 
2004). 

In the PID control there are 3 gain pieces that work on errors, namely Kp, Ki, and Kd. 
These three values greatly influence the optimal level of PID control (Priyambodo et al., 
2004). Tuning Kp, Ki, and Kd can be done manually to achieve the appropriate system 
response. The traditional method that is popularly used in PID is Ziegler Nichol. The 
weakness of the Ziegler Nichol method is that it often does not provide optimal tuning. 
Furthermore, this method easily generates high overshoot waves. Besides using the Ziegler 
Nichol method, there are also methods that use the intelligence approach. This approach uses 
optimization methods based on natural evolution. 

Evolutionary computation is a method that from the very beginning of its creation 
mimicked Charles Darwin's theory of evolution, "Survival of the fittest" (Priyambodo et al., 
2004). One method that applies optimization techniques based on natural evolution is through 
the process of mutation, crossovers and selection is Genetic algorithm (Sivanandam and 
Deepa, 2008). The unknown parameter value for PID is definitely coded as a chromosome in 
binary code. The binary code results are then evaluated using objective functions based on 
costs in time domain (Jones and Oliveira, 1995). There are also several studies that implement 
genetic algorithm methods in quadrotor PID control systems (Abdollahi et al., 2015; Bolandi 
et al., 2013; Li et al., 2014). But the disadvantage of applying the Genetic Algorithm rule is 
that it can reject important information contained in other individuals and can cause premature 
convergence. 

Premature convergence is a situation where Genetic Algorithm often cannot get a global 
optimal solution but is trapped into a local optimal solution. This situation is also supported 
by the application of static probability values on the crossover rate and mutation rate. Based 
on the problem, this study tries to apply different rules for determining the crossover rate and 
mutation rate to improve the search process for the values of Kp, Ki, and Kd which are 
optimal for quadrotor transfer functions. Determination of crossover rate and mutation rate is 
adjusted to the fitness value of each individual. This condition can cause the probability of 
crossover rate and mutation rate to decrease according to increasing generation. The decrease 
in crossover rate and mutation rate is applied when it is found that the individual fitness value 
is higher than the average individual fitness in previous evolution. Crossover rate and 
mutation rate will return static value if it is found that individual fitness values are lower or 
equal to the average individual fitness in previous evolution. 
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METHODOLOGY 

The method used in tuning PID in this study is to use the Genetic Algorithm that has been 
modified so that it is not trapped into a local optimal solution or called Fast Genetic 
Algorithm. To be able to apply this method, then do the various initial stages to get the data 
needed, namely: 

Modeling Analysis 

Newton-Euler approach to determine the dynamics of quadrotor, namely through equations 
(1) (2) (3) (4): 


^ = V 

(1) 

mv = / 

(2) 

R = Rti 

(3) 

lil = —fix IQ. + T 
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The quadrotor inertia moment is described in the equation (5)(6)(7): 
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Where M is the quadrotor middle box mass, m is the brushless motor mass, R is the 
distance between the ends of the middle part of the quadrotor to the center of mass, r is the 
brushless motor radius, h is the brushless motor height, H is the height of the quadrotor 
center, and 1 is the length between the ends of the quadrotor to the center of mass. The final 
equation for angular acceleration (j)", 0", and \|/"is linearized based on the drift state into the 
equation (8)(9)(10): 
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Transfer function 

The block diagram of Laplace's domain transfer function pitch angle, roll angle, and yaw 
angle in the quadrotor is in accordance with Figure 1. 
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Figure 1. The PID Controller Inserted Block Diagram: (a) Roll Angle, (b) Pitch Angle, (c) Yaw Angle 


Fast Genetic Simulation Algorithm 

In Fast Genetic Algorithm, the crossover (Pc) probability value decreases with increasing 
generation. This causes the population size in each generation to decline. But it can speed up 
the computing process. Equation (11) describes the determination of the probability value 
crossover (Pc). 


Pc= P,i-((m-l)*P,i/M) (11) 

Where: 

Pci = The probability of a specified crossover is 1 

m = Generation evolution 
M = Total Generation 

The mutation process that will be used in this study is a mutation for binary numbers. A 
simple way to do binary mutations is to replace one or more gene values from chromosomes. 
The mutation method is known as Flip Bit. The mutation process also depends on the 

probability of mutation (Pm). In Fast Genetic Algorithm, the probability value of mutation 

(Pm) decreases every generation. The goal is also to speed up the computing process. 
Equation (12) describes the determination of the mutation probability value (Pm): 


Pm = Pmi - ((rn - 1) * P^i/M) (12) 

Where: 

Pmi = The mutation probability is set to 0.01 

m = Generation evolution 

M = Total Generation 

Overall, the PID optimization simulation using Fast Genetic Algorithm is shown in Figure 2. 
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Figure 2. Simulation of PID Optimization with Fast Genetic Algorithm 

In addition, other useful methods were provided (Istanto and Manggau, 2018; Lamalewa 
and Maulany, 2018; Latuheru and Sahupala, 2018; Waremra and Bahri, 2018). 

RESULTS 

Simulation Implementation 

The difference between ordinary Genetic Algorithm and Fast Genetic Algorithm lies in the 
application of crossover and mutation probabilities. In Genetic Algorithm, the probability of 
crossover and mutation is set statically. In Fast Genetic Algorithm, the crossover and 
mutation probability is adjusted differently for each generation. The simulation code is 
described in Figure 3 and Figure 4. 
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%CROSSOVER 

1 Nmate=size(mpoll); 

2 random=rand(1,Nmate(1)); 

3 parentlnd=[]; 

4 j=l; 

5 for i=l:Nmate(1) 

6 if random(i) < co_rate(generasi) 

7 parentind (j)=i; 

8 j=j+l; 

9 end; 

10 end 

11 % pairs chromosomes and performs crossover 

12 for ic=l:2:length(parentind) 

13 if ic>=length(parentind) 

14 break; 

15 else 

16 Nmate=size (mpoll); 

17 cross=ceil((L-1)^rand()); 

18 mpoll(Nmate(1)+1,1:cross)=mpoll(parentind (ic),1:cross); 

19 mpoll(Nmate(l)+l,cross + l:L)=mpoII(parentind (ic+1),cross + 1:L); 

20 mpoll(Nmate(1)+2,1:cross)=mpoII(parentind (ic+1),1:cross); 

21 mpoll(Nmate(1)+2,cross + 1:L)= poll(parentind (ic),cross + 1:L); 

22 end 

23 end 


Figure 3. Code of the Fast Genetic Algorithm Crossover Program 


%MUTATION 

1 mutrate=0.001+([0:1:n_iter]^(0.01/n_iter)); 

2 Nmuts=ceil (mutrate"^ (Npop-1) "^L) ; 

3 Nmate=size(mpoll); 

4 for ic=l:Nmuts(generasi) 

5 ix=ceil (Nmate (1) "^rand) ; 

6 iy=ceil(f^rand); 

7 mpoll(ix,iy)=l-mpoll(ix,iy); 

8 end%ic 


Figure 4. Mutation Program Code Fast Genetic Algorithm Simulation Test 

Simulation Test 

When Kp is 0.0869, Ki is 0.001 and Kd of 0.0196 has a maximum peak value of 2.9 degrees. 
This means that the quadrotor will be stable to float while still being able to keep the tilt angle 
no more than 2.9 degrees. The value of settling time obtained in these conditions is 2.59 
seconds. This means that the angle of slope will decrease until it reaches stationary state with 
a travel time of around 2.59 seconds. The number of generations to achieve this condition is 
around 34 generations. The simulation execution time for this condition is 84.7093 seconds. 
Impulse response pitch and roll angles along with progress tuning for each generation are 
shown in Figure 5. 
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Figure 5. Results of pitch and roll angle simulations with Fast Genetic Algorithm 

Comparative data on the results of aehieving minimum eost in tuning Genetic Algorithm 
and Fast Genetic Algorithm simulations are shown in Table 1. 


Table 1. Comparison of search results minimum cost 


Algorithm 

Control constant 

Cost 

Achievement 
of optimum 
generation 

Execution 

time 

(seconds) 

Kp 

Ki 

Kd 

FGA 

0.0869 

0.001 

0.0196 

5.4941 

34 

84.7093 

SGA 

0.0948 

0.0126 

0.0289 

5.4938 

99 

150.068 


In Table 1 it is known that there is almost no difference between the minimum cost 
generated through the PID tuning simulation using both the Genetic Algorithm and Fast 
Genetic Algorithm. But there is an increase in the achievement of optimum generation by Fast 
Genetic Algorithm, which is the generation when the minimum cost is found: 

/34 - 99\ 

abs f ——j X 100% = 65.66% 

There was also an increase in simulation execution time between tuning simulations using 
the Standard Genetic Algorithm and Fast Genetic Algorithm, where the use of the Fast 
Genetic Algorithm method was faster than: 


abs 


r84,7093 - 150,068\ 
^ 150,068 / 


X 100% = 43.55% 


Dynamic Testing 

The test results are carried out by flying quadrotor at a certain height with a pwm value of 
more than 1500. The PID constant used is the PID constant tuning using the Fast Genetic 
Algorithm method. The response you want to know is the peak value and settling time. 

The pitch and roll angle response uses the constants Kp, Ki and Kd obtained from the 
simulation using Fast Genetic Algorithm shown in Figure 6 and Figure 7. 
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Figure 6 Response of dynamic test pitch angle with PID results of tuning Fast Genetic Algorithm 



Time (KlOOms) 


Figure 7. Dynamic roll angle test response with PID as a result of tuning Fast Genetic Algorithm 


In Figure 6 and Figure 7, the desired pitch and roll angle is 0 The peak value of the 
dynamic test results indicate that the quadrotor has been able to fly by keeping the pitch and 
roll angle of less than 4 degrees. These results differ slightly from the results of the simulation 
response which state the highest slope angle of 2.9 degrees. 

CONCLUSIONS 

The conclusions that can be drawn from the results of this study are as follows: 

• The results of comparison of PID tuning simulation using fast Genetic Algorithm compared to 
using the Standard Genetic Algorithm at pitch and roll angles, shows that Fast Genetic 
Algorithm has increased attainment of a minimum cost faster 65.66%. While Fast Genetic 
Algorithm simulation execution time increased 43.55% faster than the Standard Genetic 
Algorithm. 

• The peak value or slope of the quadrotor angle obtained from the dynamic test results with the 
simulation experience different values, but relatively small. 
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